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Abstract 

The dielectric constant of the incipient ferroelectric EuTi03 exhibits a sharp 
decrease at about 5.5K, at which temperature antiferromagnetic ordering of the 
Eu spins simultaneously appears, indicating coupling between the magnetism and 
dielectric properties. This may be attributed to the modification of the soft-phonon 
mode, Tiu, which is the main contribution to the large dielectric constant, by the 
Eu spins(7/i£ per Eu). By adding the coupling term between the magnetic and elec- 
trical subsystems as —g Y Yl if^i • Sj we show that the variation of the frequency 

I (iJ) 

of soft-phonon mode depends on the spin correlation between the nearest neighbors 
Eu spins and is substantially changed under a magnetic field. 
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I. INTRODUCTION 

The coupling between magnetism and transport properties in metals produces vari- 
ous interesting phenomena such as colossal magnetoresistance in perovskite manganites. 
Hund's rules dominate in this class of compounds, where the localized spins in the t 2g 
state strongly couple with the itinerant electrons in the e g state. If a magnetic field is 
applied under such condition, the configuration of the localized spins changes, giving rise 
to a large magnetoresistance. Correspondingly, it may be desirable that coupling between 
magnetism and dielectric properties exist in certain magnetic insulators where the con- 
figuration of the localized spins to some extent affects the dielectric properties. A good 
example is yttrium magnanite YMn0 3 . An inverse s-shaped anomaly in both dielectric 
constant and loss tangent has been detected near its magnetic ordering temperature, due 
to the frustration on the triangular lattice. 1 ' 2 ' 3 EuTi0 3 is another good example. Magnetic 
measurements show that EuTi0 3 is one of the few antiferromagnets with a positive Curie- 
Weiss temperature. 4 It has been reported 5 that the dielectric constant exhibits a sharp 
decrease at its Neel temperature 5.5K. The crystal structure of EuTi0 3 is presented in 
Fig.l, which comes from the results presented in Ref.5. It exhibits a typical structure of 
perovskite titanates, where Eu ions are responsible for the antiferromagnetic ordering and 
its spin vector lies in the ah plane while the atomic motion of the soft-phonon mode with 
the T lu symmetry is along the c axis. 5 From the neutron-diffraction results EuTi0 3 has 
a type-G antiferromagnetic structure i.e., the six nearest-neighbor Eu ions to a given Eu 
ion have opposite spins while the twelve next-nearest-neighor Eu ions have parallel spins. 
In contrast to most of the perovskite titanates RTi03(R=rare earth), which usually bear 
trivalent R and Ti, EuTi0 3 has a magnetic divalent Eu because of the stable 4/ 7 electron 
configuration and accordingly a tetravalent Ti (3d ). 

X-ray absorption near edge structure 6 has shown the absence of lattice distortion from 
the ideal perovskite structure in EuTi0 3 as well as in SrTi0 3 , which is an incipient ferro- 
electric. EuTi0 3 is very similar to SrTi0 3 as an incipient ferroelectric, demonstrating no 
frequency dispersion, which is verified through dielectric measurements. 5 Incipient ferro- 
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electrics are quite different from normal ferroelectrics as far as their dielectric properties 
are concerned. The dielectric constant increases with decreasing temperature, gradually 
deviating from the Curie- Weiss law, and saturates at a very low temperature. From the 
soft mode point of view, a normal incipient ferroelectric has a polar soft mode, but this 
polar soft mode does not freeze even in the neighbourhood of OK because of the stabiliza- 
tion of the paraelectric phase by the quantum fluctuations. 7 EuTiOs, however, behaves 
differently at low temperatures. Its dielectric constant shows a sharp decrease at its Neel 
temperature 5.5K. Furthermore the experimental dielectric data of EuTiOs shows good 
agreement between the behavior of the dielectric constant and the calculated spin corre- 
lation (^Si ■ between nearest neighbors Eu spins under different magnetic fields. It is 
enlightening to consider the possible relations between the magnetic and electrical sub- 
systems and we should therefore take the coupling interaction between the magnetic and 
electrical subsystem into account in such systems at low temperatures. As a matter of 
fact, for EuTiOs the superexchange between the nearest neighbours takes place through 
the mediation of the O 2p state. Moreover the O ion plays a role in the soft-phonon mode 
with Tip symmetry in the center of the Brillouin zone. Therefore the variation of the 
Eu spin configuration may affect the behaviour of the soft-phonon mode. Furthermore 
it is clear that the soft-phonon mode with low frequency is responsible for the dielectric 

properties via the empirical formula Xq ~ — m EuTiOs. Therefore a plausible 

m(ui — uj 1 ) 

speculation is that the variation of the spin configuration at Eu site near and below T N 
will affect the behavior of the soft-phonon mode and lead to special dielectric properties. 

X. S. Gao et al. 8 utilized Monte Carlo simulations on the basis of the Ising-DIFFOUR 
model 9 to investigate the phase transition in the two dimensional ferroelectromagnetic 
lattice where the spin moment and polarization interact. For EuTi0 3 , however, things 
are different. First EuTi0 3 is a typical example where the Heisenberg model applies in- 
stead of the Ising model, since the magnetic ions are in a s-type orbital and the exchange 
interaction is isotropic. Secondly, based on available experimental data for EuTiOs, Kat- 
sufuji and Takagi 5 proposed that the dielectric constant and spin correlation (^Si ■ 
are related as e = eo ^1 + a (^Si ■ Sj^i which differs from the coupling term previously 
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suggested in Ref.8. Thus it is reasonable to add a term in the Hamiltonian of the form 
g^2 ^ qf Si ■ Sj to describe the coupling interaction between the electrical and mag- 
netic subsystems. On such an assumption, we obtain analytical relations between the 
soft-phonon mode and the spin correlation as well as between the dielectric constant and 
the spin correlation, which are in good agreement with the experimental results. Our 
results show that the squared frequency of the soft-phonon mode is related to the spin 
correlation (^Si ■ S^. We also discuss the behaviour of the soft-phonon mode at different 
temperatures and under different magnetic fields. 

In our present paper we first deal with the electrical subsystem, including electrical- 
magnetic coupling, using the perturbation theory. In order to explore the soft-phonon 
behavior, we apply the soft mode theory under the mean-field approximation. Secondly, 
to further investigate the temperature and magnetic field dependence of the frequency of 
the soft-phonon mode, we study the magnetic spin correlation, which influences the behav- 
ior of the soft-phonon mode. Finally we give our numerical results and detailed discussion. 

II. MODELS AND ANALYSIS 

We consider the Hamiltonian for a three dimensional cubic system with periodic bound- 
ary conditions. Two parameters, qi and Si are proposed here to represent the electrical 
displacement at Ti site and the Heisenberg spin at Eu site, respectively. In addition, the 
coupling interaction between the electrical and magnetic subsystems is taken into account. 
Following Janssen's model 9 the Hamiltonian for this system can be presented as following: 

H = H e + H m + H me , (1) 

where H e denotes the Hamiltonian for the electrical subsystem, H rn is the Hamiltonian for 
the magnetic subsystem and H me the coupling interaction between the two subsystems. 
Here the soft mode theory under the mean-field approximation is applied to solve H e . In 
treating H e , the Eu and O ions are fixed and only Ti ions vibrate as was supposed by 
Slater. 10 According to the soft mode theory, the atoms vibrate like a hamonic oscillator 
with small additional anharmonic part. That is to say q t at each site is subjected to a 
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single-well potential. Thus H e in the absence of electric field can be written as following: 

^ e = Eill + ^)M££«> (2) 



2 

i v 



where pi is the momentum at site I, V (qi) is the localized potential function and - vwqiqv 

i v 

denotes the two-body interaction potential. We use the simplest form to describe the an- 
harmonic motion at Ti site. Thus 

V(qi) = \mQlqf + ^qf, (3) 

1 1 
where flo is the inherent frequency, —mVt^q^ is the harmonic part and -jqf the anharmonic 

one. 

The interaction between Eu 4/ spins originates from (i)superexchange(antiferromagnetic) 

through the O 2p state and also from (ii)indirect exchange through the Eu 5d state. So 
7 

H m with S — - can be written as : 
2 

H m = J2 J i& ■ % + £« • S 3 -Y^h-Si, (4) 

where Si is the Heisenberg spin at Eu site i, and the exchange integrals J\ = 0.037A; B K 
and J2 = — 0.069A;b K to produce = 5.5K represent the antiferromagnetic coupling 
between the nearest neighbours and the ferromagnetic coupling between the next nearest 
neighbours, respectively. Here h is the external magnetic field parallel to the z axis. 
and denote the once summation over the next nearest neighbors and the nearest 
neighbors, respectively. 

As we pointed out in the previous section we propose the coupling term between the 
magnetic and electrical subsystems as: 

« <^> 

where g is the coupling coefficient. Here only the quadratic term of qi is taken because of 
the time reversal and space reversal symmetries. Here I represents the lattice site of Ti 
ion, % and j denote distinct nearest neighbor Eu ions to a given Ti ion. 



First we deal with the electrical subsystem, considering electrical-magnetic coupling. 
Under the mean-field approximation, the single-ion Hamiltonian can be written as: 



H? = Hf + Hr = ^~ 
1 11 2m 



\-mnl - z x g (sl ■ S*j^)q? + --yqf - ^ v w (q v ) q h 



(6) 



where Z\ = 12 is the number of the spin correlation that will directly affect the motion of 
Ti in one unit cell and ^7<?f can be regarded as the small perturbation part. To facilitate 

— -2 ( Si ■ Sj ) and flo(g) = ^o\/l — b. The magnetic 



the further discussion, we define b 



effect on the electrical subsystem is embodied by fio(flO- Thus 

H * = 2m~ + \ mn ^ q i + \ iq ' - E ^' «• 



(7) 



The frequency of the soft-phonon mode is given through soft-mode theory: 



.li 



where 



Vn 



E 



exp 



-iq-{R x -R v ) , 



(8) 



(9) 



is the interaction potential at the wave vector q. The key point is how to calculate the 
fluctuation of displacement ot = (qi — (qi)) 2 - In incipient ferroelectrics, quantum effects 
are very important at very low temperatures, so quantum statistics must be adopted. At 
the same time we notice that no ferroelectric phase transition takes place in the incipient 
ferroelectrics and the average displacement is zero. Thus we use quantum statistics and 
perturbation theory to obtain the total displacement fluctuation at very low temperatures, 
which contains both quantum fluctuations and thermal fluctuations: 



°t = (qi ~ (qi)Y 



-coth{— — —J — 
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COth { ^T ) + 1 



(10) 



2mn (g)~""" x 2k B T ' 16m 3 ^(#) 
For EuTiC>3 the soft-phonon mode is in the center of Brillouin zone((f = 0), therefore 
according to the single mode theory under the mean-field approximation the frequency of 
the soft-phonon mode can be written: 

37 



muj ?=0 (g) = mSl Q (g) + 3<ya T - u$ =Q = 



mQl(g) 
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where 



k B To = ^^-[u^-m^(g)]. 



(12) 



Obviously if the coupling effect is not included, i.e., g — 0, Eq. (11) is just the expression of 
the frequency of the soft-phonon mode for normal incipient ferroelectrics. The frequency 
of the soft-phonon mode does not go to zero even close to OK. In order to directly exhibit 
how the spin correlation has effect on the behavior of the soft-phonon mode we arrange 
Eq. (11) as following: 



mwl Q (g) « 77iwLo(0)(l - a I Si ■ S. 



(13) 



where 



a 



24g 



o2 ( m Q 

m\l n —coin — — — 

4mQ \2k B T 



+ 



2kBTsinh 2 



(14) 



2k B T / 

Here ^=0(0) is the frequency of the soft-phonon mode as the coupling interaction is not 
considered and it is temperature independent at very low temperatures. Because the cou- 
pling interaction only exists near the Neel ordering temperature, a can be approximately 
written as 



a 



24g 



mW^ (0) 



mill 



4mf2 



The static susceptibility, which embodies the dielectric constant, is given by: 

1 



X 



l + a(S i .S j ))= X (0)(.l + <x{Si'S j )) 



(15) 



(16) 



mul =0 (0) 

where x(0) is the static susceptibility when the coupling interaction is not included. Thus 



the dielectric constant can be derived from obove: 



e = e(p)(l + P(S i -S j )) 



(17) 



where 



e(0)=A + X (0), 



(18) 
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and 

0= X( ° )a (19) 

p A+ X (oy [ } 

Here A is the background dielectric constant. Eq. (17) is in good agreement with the 
dielectric properties fitted from the experimental data 5 . 

It is easy to show that there exists immediate relations between the behaviour of the 
soft-phonon mode and spin configuration. To investigate how the frequency of the soft- 
phonon mode varies under different conditions of spin correlation, we take into account the 
spin correlation in the coupling system. Then Hamiltonian for the magnetic subsystem 
can be written: 

H M = H m + H me = J2(Ji-9j2li^-^ + J ^J2^-^-J2 J '-^- ( 2 °) 

This magnetic subsystem can be divided into two sublattices a and b. Thus 

H M = H a + H b . (21) 

where H a and H b are the Hamiltonians belong to a and b sublattices, respectively. Under 
the mean-field approximation 

Ha = H* a J2 SX ai + HZ aJ2 S ai> ( 22 ) 



and 



where 



H b = H^Su + H z b J2S z bl , (23) 



H% b) = 6(Ji - z 2 g (qf)) (£% a) ) + 12J 2 (£% b) ) , (24) 
H z a{b ) = ~h + 6( Ji - z 2 g (qf)) (Sfa) + 12J 2 (S z a(b) ) . (25) 

Here x, z denote the spin components and z 2 = 4 is the number of Ti that will directly 
affect the spin configuration at each pair of correlation. Thus the spin components of the 
different sublattices can be obtained: 

™ E (2t - l)*nh&f±y/{H S ) 2 + (Htf) 

( S «) = 1 = v —> ( 26 ) 

+ m 2 E cosh(^y/(Htf + {Htf) 



i=i 



(St) 



(S x a ) 




(27) 



(28) 



(29) 



The spin correlation between the nearest neighbours can be approximately decoupled as: 



For a complex system where coupling exists between the magnetic and electrical subsys- 
tems one should consider the subsystem results obtained by including the coupling term 



III. RESULTS AND DISCUSSION 

Using Eq. (17)-Eq. (19) and Eq. (22)-Eq. (30) we fit the experimental results 5 in 
the absence of magnetic field. We not only get the parameter To = — 25K, which is the 
same as that in Ref.5, but also we get the coupling coefficient g as 9.42 x 10~ 7 , which is 
doubtlessly small but plays an important role on the behavior of the soft-phonon mode. 
The parameters f2 and 7 are also calculated as: f2 ~ 2.13 x 10 13 and 7 3.59 x 10 19 . 
Now we use these parameters to study the effect of coupling on the frequency of the 
soft-phonon mode. 



the presence of electrical-magnetic coupling is shown in Fig. 2. Variation of the spin 
configuration can be realized by varying the temperature and the external magnetic 
field. The spin correlation does not change significantly if the coupling interaction is 
not neglected(g — > 0). That is to say, the electrical effect is not strong enough to sub- 
stantially reduce the spin correlation between the nearest neighbours. 

But the magnetic influence on the behavior of the soft-phonon mode is rather evident. 
The temperature dependence of the frequency of the soft-phonon mode under different 




(30) 



l (ij) 




versus temperature under different magnetic field in 
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magnetic field is presented in Fig.3. We focus on the soft mode properties for T < 30K, 
where the frequency of the soft-pnonon mode as well as the dielectric constant shows little 
temperature variation if the electrical-magnetic coupling is neglected(<7 = 0; solid line). 
When the electrical-magntic coupling is included, the soft-phonon frequency increases 
remarkably below the antiferromagnetic ordering temperature (about 5.5K for h — 0), 
instead of the soft mode behavior of normal incipient ferroelectrics. One can infer that 
the antiferromagnetic ordering stiffens the soft phonon mode. Consequently, the dielectric 
constant decreases below the Neel temperature. However, when higher magnetic fields 
are applied, sufficient to induce ferromagnetic ordering, softening of the soft-phonon mode 
can be observed. Such softening is responsible for the increase of the dielectric constant 5 . 
To summarize, our results suggest that the soft-phonon mode is modified by the Eu spins. 
The configuration of Eu spins determines the hybridization between the Eu orbitals and 
the O 2p orbitals and then gives rise to the modification of the behaviour of the soft- 
phonon mode that contains the Ti - O stretching motion. 

The soft-phonon frequency versus applied magnetic field was also computed(Fig.4). 
We can see that for T < 20K, the frequency of the soft-phonon mode decreases with 
increasing magnetic field, but the exact behavior of u(h) depends on T. At 2K, u(h) 
decreases sharply as h increases for h < 1.5T. For h > 1.5T, u(h) decreases more slowly 
as h increases, owing to lesser variation of the spin correlation (^Si ■ S^. For T = 4K, 
ui{h) shows less variation, however a kink in u>(h) still appears. At 5.5K the curve is more 
smooth and only a very slow decrease with the increasing magnetic field is observed. 
When a high magnetic field is applied, the spin correlation under the ferromagnetic or- 
dering lowers the frequency. The lower the temperature the stronger the effect. For low 
magnetic field, however, the antiferromagnetic ordering dominates and the frequency of 
soft-phonon mode is enhanced relative to u(0). The lower the temperature the greater the 
enhancement. For T > 20K, at which the spin correlation becomes minuscule, u(h,T) 
is nearly independent of h and increases as T increases. This is the typical behavior for 
normal incipient ferroelectrics within the framework of classical statistics. At such high 
temperatures, the spin correlation is too weak to significantly affect the soft mode. Under 
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such conditions the frequency of the soft-phonon mode in quantum statistics reduces to 
the classical Curie- Weiss result: cu oc (T — T )^ 



IV. SUMMARY: 

EuTi0 3 is regarded as both a typical type-G antiferromagnet with S — \ at Eu 2+ site 
and an incipient ferroelectric with the perovskite crystal structure. This special property 
makes the frequency of the soft-phonon mode as well as the dielectric constant deviate 
from the behavior in nomal incipient ferroelectric, resulting from coupling between the 
magnetic and electrical subsystems. From our discussion above, variation of the spin 
correlation l^Si • S^j causes a remarkable variation of the soft-phonon frequency, which 
is regarded as the main reason for the special dielectric properties of EuTi0 3 . 
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CAPTION OF FIGURES 



Fig.l: The crystal structure of EuTi0 3 . White circles and a black circle represents O 
and Ti, respectively. The white thick arrows are the antiferromagnetic ordering of Eu 
spins. The black thin arrows exhibit the atomic motion of the soft-phonon mode. 



Fig. 2: The temperature dependence of pair correlation (^Si ■ Sjj between the nearest 
neighbours at Eu 2+ . (g = 9.42 x KT 7 ) 

Fig. 3: The temperature dependence of the frequency of soft-phonon mode under dif- 
ferent magnetic fields, (g = 9.42 x 10~ 7 ) 

Fig4: The frequency of the soft-phonon mode versus magnetic field at different tem- 
peratures, (g = 9.42 x 10 -7 ) 
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